function uswm_void_fraction_analysis()

       close all
       clear all
    us=0;
    wm=0;
    filename = ['.\Results_alpha\Workspace_analysis_timeseries_hist_PSD__selection_void_US1.mat'];    
    if (exist(filename,'file'))          %If there were treatments of LVM files_list        
        filevariables = {'matrix_mean','alpha_gregory_scott','alpha_mattar_gregory'};             
        load(filename, filevariables{:});        
        matrix_us=matrix_mean;    
        clear matrix_mean;
        us=1;
    end
    
    filename = ['.\Results_alpha\Workspace_analysis_timeseries_hist_PSD__selection_void_WM1.mat'];    
    if (exist(filename,'file'))          %If there were treatments of LVM files_list        
        load(filename,'matrix_mean');
        matrix_wm=matrix_mean;
        wm=1;
    end
    if(us==1 && wm==1)
    
        figure(5);                                    %print graph mean_void_fraction vs JG/J
        whitebg('w')    %White background
        set(gcf, 'color', 'w');
        set(gca, 'color', 'w');
        ftype = fittype('a+b*x');
        fitus = fit(matrix_us(:,10), matrix_us(:,1),ftype);
        fitwm = fit(matrix_us(:,10), matrix_wm(:,1),ftype);
        
        
        plot(fitus,'k--',matrix_us(:,10), matrix_us(:,1),'ko'); %void fraction vs JG/J        
        %plot(fitus,'k--',matrix_us(:,10), matrix_us(:,1))        
        hold on        
        %plot(matrix_us(:,10), matrix_us(:,1),'ko','MarkerFaceColor','k'); %void fraction vs JG/J        
        plot(fitwm,'b--', matrix_us(:,10), matrix_wm(:,1),'bo'); %void fraction vs JG/J
        %Plots Gregory and Scott/Mattar and Gregory graphic
        plot(matrix_us(:,10), alpha_gregory_scott,'ks','MarkerEdgeColor','r'); %void fraction vs JG/J
        plot(matrix_us(:,10), alpha_mattar_gregory,'k^','MarkerEdgeColor','g'); %void fraction vs JG/J
        grid on
        %plot(matrix_us(:,10), matrix_wm(:,1),'k^','MarkerFaceColor','b'); %void fraction vs JG/J
        legend('Measurement Data US','Fitted curve US','Measurement Data WM','Fitted curve WM','Gregory and Scott','Mattar and Gregory','Location', 'Best');
        xlabel('JG/J'); 
        y_label_aux = ['Mean \alpha'];
        ylabel(y_label_aux);
        file_path_output = ['.\Results_alpha\Figure_alpha_J_models_USWM.fig'];
        saveas(gcf,file_path_output); 
        file_path_output = ['.\Results_alpha\Figure_alpha_J_models_USWM.tif'];
        saveas(gcf,file_path_output); 
        
        figure(10);                                    %print graph mean_void_fraction vs JG/J
        whitebg('w')    %White background
        set(gcf, 'color', 'w');
        set(gca, 'color', 'w');
        ftype = fittype('a+b*x');
        fitus = fit(matrix_us(:,10), matrix_us(:,1),ftype);
        fitwm = fit(matrix_us(:,10), matrix_wm(:,1),ftype);
        
        
        plot(fitus,'k--',matrix_us(:,10), matrix_us(:,1),'ko'); %void fraction vs JG/J        
        %plot(fitus,'k--',matrix_us(:,10), matrix_us(:,1))        
        hold on        
        %plot(matrix_us(:,10), matrix_us(:,1),'ko','MarkerFaceColor','k'); %void fraction vs JG/J        
        plot(fitwm,'b--', matrix_us(:,10), matrix_wm(:,1),'bo'); %void fraction vs JG/J        
        grid on
        legend('Measurement Data US','Fitted curve US','Measurement Data WM','Fitted curve WM','Location', 'Best');
        xlabel('JG/J'); 
        y_label_aux = ['Mean \alpha'];
        ylabel(y_label_aux);
        file_path_output = ['.\Results_alpha\Figure_alpha_J_USWM.fig'];
        saveas(gcf,file_path_output); 
        file_path_output = ['.\Results_alpha\Figure_alpha_J_USWM.tif'];
        saveas(gcf,file_path_output); 
%         %%
%          figure(11);                                    %print graph mean_void_fraction vs JG/J
%         whitebg('w')    %White background
%         set(gcf, 'color', 'w');
%         set(gca, 'color', 'w');
%        
%         plot3(matrix_us(:,7),matrix_us(:,8),matrix_us(:,1),'o','MarkerFaceColor','w')
%         %plot(fitus,'k--',matrix_us(:,10), matrix_us(:,1))        
%         hold on        
%         
%         plot3(matrix_us(:,7),matrix_us(:,8),matrix_wm(:,1),'ro','MarkerFaceColor','w')
%         %plot(matrix_us(:,10), matrix_us(:,1),'ko','MarkerFaceColor','k'); %void fraction vs JG/J        
%         
%         grid on
%         %legend('Measurement Data US','Fitted curve US','Measurement Data WM','Fitted curve WM','Location', 'Best');
%         %xlabel('JG/J'); 
%         %y_label_aux = ['Mean \alpha'];
%         %ylabel(y_label_aux);
%         %file_path_output = ['.\Results_alpha\Figure_alpha_J_USWM.fig'];
%         %saveas(gcf,file_path_output); 
%         %file_path_output = ['.\Results_alpha\Figure_alpha_J_USWM.tif'];
%         %saveas(gcf,file_path_output); 
       %%     
        figure(7)
        plot_referenceXexperimental(alpha_gregory_scott',matrix_us(:,1))
        grid on
        xlabel('gregory scott \alpha'); ylabel('ultrasound \alpha');    %Lambels the axis
        file_path_output = ['.\Results_alpha\Figure_alpha_J_gregscott_US.fig'];
        saveas(gcf,file_path_output); 
        file_path_output = ['.\Results_alpha\Figure_alpha_J_gregscott_US.tif'];
        saveas(gcf,file_path_output); 
        
        figure(8)
        plot_referenceXexperimental(alpha_mattar_gregory',matrix_us(:,1))
        grid on
        xlabel('mattar gregory \alpha'); ylabel('ultrasound \alpha');    %Lambels the axis
        file_path_output = ['.\Results_alpha\Figure_alpha_J_mattargregory_US.fig'];
        saveas(gcf,file_path_output); 
        file_path_output = ['.\Results_alpha\Figure_alpha_J_mattargregory_US.tif'];
        saveas(gcf,file_path_output); 
        
        
        figure(9)
        plot_referenceXexperimental(matrix_wm(:,1),matrix_us(:,1))
        grid on
        xlabel('wire-mesh \alpha'); ylabel('ultrasound \alpha');    %Lambels the axis
        file_path_output = ['.\Results_alpha\Figure_alpha_J_USWM.fig'];
        saveas(gcf,file_path_output); 
        file_path_output = ['.\Results_alpha\Figure_alpha_J_USWM.tif'];
        saveas(gcf,file_path_output); 
        
    end
    
    

return